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Abstract 

We study J ftp suppression in Pb+Pb collisions at CERN-SPS energies in hadronic 
matter with energy- and temperature-dependent charmonium dissociation cross sec- 
tions calculated in the quark-interchange model of Barnes and Swanson. We find 
that the variation of J/ip survival probability from peripheral to central collisions 
can be explained as induced by hadronic matter absorption in central collisions. 
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I. INTRODUCTION 



Following the suggestion of Matsui and Satz [1] to probe quark-gluon plasma (QGP) 
with heavy quarkonium, the NA50 Collaboration investigated J/ip suppression in high- 
energy heavy- ion collisions. Anomalous J /if) suppression has been observed in Pb+Pb 
collisions at 158 A GeV [2-3]. It is of interest to investigate whether deconfined matter 
is created in these collisions. Many different mechanisms have been proposed to explain 
the phenomenon. It has been suggested that the suppression was due to a new phase, 
the QGP [4], a change of the equation of state due to the QCD phase transition [5], the 
melting of charmonia in QGP [6], or a percolation deconfmement [7]. It has also been 
suggested that the anomalous suppression comes from the absorption by various comovers 
produced in the Pb+Pb collisions [8-14]. 

In addition to the general suppression over a large region of transverse energy E T , 
which is related to the impact parameter b, the data for J/ip suppression in the large 
transverse energy region with E T > 100 GeV are also of special interest. This region 
corresponds to the most central Pb+Pb collisions. The possible presence of a rapid drop 
of the ratio of J/ip to Drell-Yan production cross sections, B(J /ip)a(J /ip) / a(DY) , has 
been interpreted in a comover model as a result of multiplicity fluctuation, subject to 
a large uncertainty of the measured data [15]. It has also been suggested that large 
E T fluctuations lead to increasing probability for deconfmement and therefore possible 
complete suppresssion of J/ip [16]. A definitive understanding on the nature of J/ip 
suppression in Pb+Pb collisions is still lacking. 

The absorption of J / ip by comovers depends on the J / ip dissociation cross sections in 
collision with hadrons. The J / ip dissociation cross sections in collision with hadrons have 
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been considered previously in several theoretical studies, but the predicted cross sections 
show great variation at low energies, largely due to different assumptions regarding the 
dominant scattering mechanism. Kharzeev and Satz and collaborators [17] employed the 
parton model and perturbative QCD "short-distance" approach of Bhanot and Peskin 
[18] and found remarkably small low-energy cross sections for collisions of J ftp with light 
hadrons. Matinyan and Miiller [19], Haglin [20], Lin and Ko [21], Song and Lee [22], 
and Navarra and collaborators [23] recently reported results for these dissociation cross 
sections in meson exchange models. These references all use effective meson Lagrangians, 
but differ in the interaction terms included in the Lagrangian. Results depend on a 
choice of meson coupling schemes and assumed form factors. Charmonium dissociation by 
nucleons has also been considered recently using similar effective Lagrangian formulations 
[24]. 

We favor the use of the known quark forces to obtain the underlying scattering 
amplitudes from explicit nonrelativistic quark model wavefunctions of the initial and final 
mesons [25, 26]. In this approach we first look for the inter-quark interaction by using 
the meson spectrum. The wave function and the interaction are then used to evaluate the 
cross section for the quark-interchange process [25, 26]. The calculated cross sections are 
first compared with experimental data of / = 2 tut phase shifts, and excellent agreement 
between theory and experimental data was obtained [25, 26]. Such an approach allows 
us to calculate all charmonium dissociation cross sections by collisions with hadrons. 
Martins, Blaschke, and Quack previously reported dissociation cross section calculations 
using essentially the same approach [27]. Two general characteristics of the cross sections 
obtained in this approach are: (1) for endothermic reaction, the cross section rises after 
the hadron energy exceeds some threshold energies and the maximum magnitude of the 
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cross section is in the mb range, (2) for exothermic reaction, the cross section is infinite 
at zero kinetic energy and decreases as the energy increases. 

The J/ip dissociation cross section depends on the meson wave functions and the 
inter-quark interaction. The latter quantities depend on temperature. By using such a 
temperature-dependent inter-quark potential inferred from lattice gauge calculations [28] , 
the temperature dependence of J/ip dissociation cross sections can be obtained. The new 
cross sections are now applied to study B(J/ip)a(J/ip) /a(DY) versus E T , to be compared 
with the CERN-SPS Pb+Pb data of the NA50 Collaboration. 

The concept of the temperature dependence of the dissociation process should be 
clarified as many different processes are involved. There are two types of thermalization 
[26,29]. First, there is the thermalization of the hadronic matter, in which the J/ip (or its 
precursor) resides. It occurs rapidly as it depends on the cross section for the scattering 
of light hadrons with light hadrons, which is of the order of tens of mb. We expect that 
the dense hadronic matter is thermalized within a time of the order of 1 fm/c after light 
hadrons are produced. (In our numerical results below, the thermalization time for the 
hadronic matter is 1.9 fm/c after the light hadrons are produced.) 

The second type of thermalization is the thermalization of the "charmonium sys- 
tem" with its surrounding medium. If a quarkonium is placed in a thermalized hadronic 
medium, there will be non-dissociative inelastic reactions between the quarkonium and 
medium particles which change a charmonium state into another charmonium state: 
h+ (QQ)jls <-> h' + (QQ)j'L'S'- These reactions lead to the thermalization of the charmo- 
nium system [26,30]. When the heavy quarkonium system is in thermal equilibrium with 
the medium, the probabilities for the occurrence of heavy quarkonium states will be dis- 
tributed according to the Bose-Einstein distribution. The time for the thermalization of 
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such a system depends on the non-dissociative inelastic cross sections. As these cross sec- 
tions are different from the pion-pion scattering cross section, the "charmonium system" 
thermalization time is different from the light-hadron matter thermalization time. 

The cross sections for non-dissociative inelastic reactions between a charmonium and 
a light hadron and their temperature dependences remain a subject of current research. 
Previous estimates of n + J/tp — > n + ip' by Fujii and Kharzeev [31] gave a cross section 
of the order of 0.01 mb. In the quark-interchange model of Barnes and Swanson [32,33], 
such a process is forbidden in the first order and takes place only in the second order. The 
higher order process presumably will lead to a smaller cross section. The non-dissociative 
inelastic cross sections need to be evaluated for the interaction of pions with different 
charmonia using different models, to provide an estimate of the time for the thermalization 
of the charmonium system with the hadronic medium. If the non-dissociative inelastic 
charmonium cross sections are of the same order as in 7r — 7r collision, then the charmonium 
system will be also thermalized and the dissociation of J/ip should be treated by the 
method of statistical mechanics. On the other hand, if the non-dissociative inelastic 
cross sections are all of the order of 0.01 mb, as estimated by Fujii and Kharzeev for 
7r + J/tp — > it + ip', the thermalization time for the charmonium system will be quite long 
and we will not need to consider a thermalized charmonium system. In the present work, 
we have not taken the process of dissociation by thermalization into account. 

For a charmonium placed in a thermal medium, the thermalized hadronic matter 
alters the interaction between the charm quark and charm antiquark and changes the 
QCD vacuum surrounding the charmonium. The threshold energy for J/ifj dissociation 
is shifted and consequently the dissociation cross section changes as a function of the 
hadronic matter temperature. The dissociation cross section is however of the order of 
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a few mb [26]. Thus, in the treatment of the collision of J/ip with the dense hadronic 
matter, the cross section is still small enough so that we can treat the n — J/ip dissociative 
collision as a two-body process. 

In Sec. II we describe the inter-quark potential used to obtain meson wave functions. 
In Sec. Ill we present the n — J/ip and 7r — Xcj dissociation cross sections and show their 
dependences on temperature. In Sec. IV charmonium suppression in Pb+Pb collisions 
is examined. Numerical calculations and results are shown and discussed in Sec. V. 
Conclusions are given in Sec. VI. 

II. TEMPERATURE-DEPENDENT POTENTIAL AND WAVE FUNCTIONS 

In a medium at high temperatures, the gluon and light quark fields fluctuate and 
the alignment of the color electric fields due to the QCD interaction is reduced by the 
thermal motion for a random orientation of the color electric fields. The "pressure" from 
the QCD vacuum to confine the quark and antiquark pair also diminishes with increasing 
temperature. The inter-quark potential is thus a sensitive function of temperature, re- 
sulting in the vanishing of the string tension at the critical temperature of deconfinement 
phase transition. 

Even at T = 0, a proper description of the heavy quarkonium state should be based 
on a screening potential, as the heavy quarkonium becomes a pair of open charm or open 
bottom mesons when r becomes very large, due to the action of dynamical quark pairs 
[34] . Recently, the central inter-quark potential has been obtained from the lattice gauge 
results of Karsch et al. [28] and analyzed by Digal et al. [35, 36] and by Wong [26]. The 
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temperature-dependent potential for T < T c can be represented by 



where a s is a running coupling constant [25]. The effective string-tension coefficient is 

b(T) = b [l-(T/T c ) 2 ]9(T c -T), (2) 

with b = 0.35 GeV 2 and the effective screening parameter 

fi(T) = fi 9(T c -T), (3) 

where /x = 0.28 GeV and the step function 9 signifies the vanishing of the string tension 
at the phase transition temperature T c = 0.175 GeV. 

The charmonium wave functions can be obtained by solving the Schrodinger equa- 
tion 

[-V • + V 12 (r,T) + A(r,T)^ JLS (f,T) = e(T)rl> JLS (r,T), (4) 

^12 

where /i^ is the reduced mass and the mass difference A(r, T) is 

A(r, T) = mi(r,T) + m 2 (r,T) - m 1 (oo,T) - m 2 (oo,T). (5) 

The energy e(T) is measured relative to two separated mesons cq and qc at r — > oo where 
c and c become open charm mesons and {mi(r, T), m 2 (r, T)} = {M C g(T), Mg g (T)}. 

The bound-state wave function and energy e have been obtained by diagonalizing 
the Hamiltonian using a nonorthogonal Gaussian basis with different widths [26]. When 
the energy e becomes positive, spontaneous dissociation occurs, subject to the selection 
rules based on the conservation of total angular momentum and parity [26]. 

III. CHARMONIUM DISSOCIATION CROSS SECTIONS 



For meson-meson scatterings A + B — > C + -D, the differential cross section is given 

by 

where t is the momentum transfer squared, s is the center-of-mass (CM.) energy squared, 
and pa is the momentum of meson A in the CM. system. The matrix element Aifi is 
calculated with the four meson wave functions of A, B, C, and D and the inter-quark 
interaction. The latter is the potential in Eq. (1) plus a spin-spin term for quark-quark 
interaction, 

A A 

v„ iQ (r) = -y • y x ' g ^ e ~*'°> (7) 

where mi(m 2 ) and Si(s*2) are individually the mass and spin of quark 1(2), and d is a 
spin-spin interaction width including relativistic corrections [26, 37]. For quark- ant iquark 

A A A A T 

potential, — • y is replaced by — — • -y. This potential generates wave functions which 
are then used to calculate the dissociation cross sections, as discussed in Ref. [26]. The 
dissociation cross section of J ftp in collision with tc varies as a function of T/T c and the 
CM. kinetic energy E KE = y/s — — m B where and m B are the masses of the 
incident mesons. We plot in Fig. 1 the ir + J ftp dissociation cross sections. Since about 
35% of J/ip comes from the radiative decay of Xcj, the dissociation cross sections of Xci 
and Xc2 in collisions with n are also important. We evaluate these cross sections which 
are plotted in Fig. 2. 

As the temperature increases, the temperature-dependent inter-quark potential be- 
comes weaker. Then the energy for separating c and c in the charmonium decreases and 
the root-mean-square radius of the meson increases. As a consequence, the dissociation 
cross section for a charmonium in collision with a pion rises with increasing temperatures. 
For various charmonia, different sizes also lead to different dissociation cross sections. At 
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T/T c PS 0.7, the peak cross section of tt + Xci or 7r + Xc2 is about 2.1 times that at T = 0. 
Since the Xc2 is less bound compared to Xci, the dissociation cross section for Xc2 is slightly 
larger. The n + Xc2 reaction becomes exothermic at a lower temperature than n + Xci- 
Comparing Figs. 1 and 2, we find that the peak value oin — J/tp dissociation cross section 
is larger than the peak values of n — Xc2 and it — Xci dissociation cross sections when T/T c 
are larger than 0.65 and 0.55, respectively. 

IV. CHARMONIUM SUPPRESSION IN Pb+Pb COLLISIONS 

Working in the nucleon-nucleon CM. frame, we assume the following scenarios to 
investigate J/ip suppression. In a nucleus-nucleus collision, hadronic matter consisting 
mainly of pions is produced at rf orm after the collision. The hadronic matter is initially 
not in thermal equilibrium. It becomes thermalized by elastic scatterings of pions at time 
r thcrm, and the expansion of hadronic matter lowers the temperature and the density until 
freeze-out at Tf z . On the other hand, cc pairs are created by the hard scattering processes 
between a projectile nucleon and a target nucleon in a very short time. The cc pair 
evolves into a charmonium or its precursor which interacts and becomes dissociated by 
colliding first with nucleons in the two colliding nuclei and subsequently with the produced 
hadronic matter. J/vf) will be dissociated by collision with pions during the time from 
Tform to Tf z . The differential cross section for J/ip production with respect to the impact 
parameter b is 

do j a, a p j o NJ /^ 

/"^thcrm 

xexp[-/ drv^a^j/^n^b, r)] 

-'''"form 
/"Tfz 

x exp[- / driv^a^j/^n^b, r)]. 

J r thcrm 
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For the production of Drell-Yan dilepton pairs, the cross section scales with the 
number of nucleon-nucleon collisions, and we have 

^M^ = a^ABT AB (b). (9) 
do 

The quantities E and p are the energy and momentum of J/ip; a™ and a^y are cross 
sections for J ftp and Drell-Yan productions in a nucleon-nucleon collision, respectively; 
0"ttJ/v> 1S the J /if) dissociation cross section in collision with n; a^j/ip is the J/tp absorption 
cross section in collision with nucleons; b A is the nucleon coordinate in the target nucleus 
A; T A and T B are nuclear thickness functions [38]; and v rel is the relative velocity of 
charmonium and ir. The thickness function for the colliding nuclei is 

T AB (b) = J db A T A (b A )T B (b - b A ). (10) 

The expression for the differential cross section for Xcj and ip' production can be obtained 
by replacing J/ip with Xcj or ip' . 

If the hadronic matter is in thermal equilibrium, the pion number multiplicities are 
obtained from the familiar momentum distribution. Before thermalization, the pion num- 
ber density has to be estimated from experimentally measured multiplicities of charged 
pions. The total pion number N n as a function of the impact parameter b for Pb+Pb 
collisions at 158 A GeV was measured by the NA49 collaboration [39]. We assume that 
the total pion number is conserved in the evolution of hadronic matter. The volume of 

— * 

hadronic matter depends on the impact parameter b as well as the proper time r, which 
is measured relative to the initial collision, 



JAR 2 A - b 2 
V(b,r) = 2r+ 1 
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x{2(R A + rW ™- \ R J tVc J - V + ^ " t}-<»> 
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Here Ra is the lead-nucleus radius, V cx ±_ is the transverse velocity for radial flow, and 7 
is the Lorentz contraction factor, 7 = y/s/2rriN, with nucleon mass mjv, and ^J~s is the 
CM. energy of a nucleon-nucleon collision. As a result of the above consideration, the 
pion number density depends on the impact parameter and proper time as 

».M = ^|. (12) 

The 7r-charmonium dissociation cross sections depend on the temperature T and the 
collision energy. Pions in thermal equilibrium have the momentum distribution f n (k) ~ 
3 e --E w /T ^^11 the pion energy E n . Integrating the whole pion momentum contribution to 
i> re ic and dividing by the pion number density, we obtain 

^- I(7 > = f tfik f nA ■» ( 13 ) 

J ( 27r )3 J '7rV ft > 

which depends on p T and Xp of charmonium and r. 

The differential cross sections with respect to E T for J/tp production is given by 

_|*./^_p fl(6 , £T) , (14) 

and the differential cross section for the Drell-Yan process is 



dE T J db 
where D(b, E T ) relates b to E T [40] 



fdPb^D(b,Er), (15) 



£>( 6> E T ) = -=L e -[ET-E T (pT/2**{b)^ (16) 

\ / 2na(b) 



with the average transverse energy 



- , N 135GeV , N 

^ = l+e»p[( t -6.7)/3] - (17) 
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and the variance 



a{b) = l + exp[(6-10)/2] - (18) 



V. RESULTS AND DISCUSSIONS 

A. Differential Cross Sections for J/ip Yield in Nucleon-Nucleon Collisions 

To generate momentum distribution of J/ip, we need to parametrize the initially 
produced J/ip distributions. The invariant differential cross section for J/tp production 
in a nucleon-nucleon collision, Ed 3 aj^/d 3 p, is factorized as 

E^^ = f(x F )g(p T ). (19) 

The parameters of /(xf) and g{pr) obtained by fitting experimental data in p — p and 
p — A collisions have been summarized in Refs. [41, 42]. We take the parametrizations 

g(p T )~(i + p- 6 , (20) 

J\ x f) ~ ; , (21) 

Xi + x 2 

where x lj2 = 0.5(^Jx 2 F + 4Mj^Js ± x F ) with the J/tp mass Mj/^. The quantity (3 
is related to the average transverse momentum (p?) by (3 = fff (pt)- The quantity 
P = 2.24 ± 0.28 GeV/c is inferred from (p T ) = 0.96 ± 0.12 GeV/c for the J/ip pro- 
duction in p-Be collisions at y/s = 16.8 GeV [43]. The parameter a = 4.95 fits the 
x F spectra of J ftp in p-Be collisions at y/s = 38.8, 31.6, and 16.8 GeV [43, 44] and 
p-Li collisions at y/s = 23.8 GeV [45]. The x F spectrum provides the average value 
(x F ) = Jo ^ dx F x f f (^x f) I Jo 6 dx F f(x F ) ~ 0.16 with 8 = Mj,^/s. The condition x ij2 < 
1 requires xf < 1-5- We use the parametrizations of /(x.f) and g(pr) to generate the 
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xp and Pt spectra of a charmonium produced in nucleon-nucleon collisions in Pb+Pb 
collisions at 158 GeV/c per nucleon. 

In Eq. (8) we assume J/ip mesons are produced uniformly in the collision region. 
The successive collisions of a projectile nucleon in the target nucleus cause the nucleon 
to lose energy gradually so that J/ip mesons are not evenly produced in space. Such an 
effect was discussed in Ref. [46]. 

B. J/ip Suppression in Hadronic Matter 

As seen in Eq. (8), the differential cross section explicitly depends on the pion 
formation time Tf orm , the pion thermalization time r t h e rm, and the hadron freeze-out time 
Tf z . We estimate the formation time Tf orm to be 1 fm/c from nucleon-nucleon collision 
data [38]. Then in the most central Pb+Pb collision, n n {b = 0,r form ) = 2.38 fm~ 3 . We 
use a freeze-out number density of ri{ z = 0.5 fm~ 3 . The initial temperature of pion matter 
is close to the QCD phase transition temperature T c = 0.175 GeV [28] and is taken as 
T thcrm = 0.16 GeV. 

In a thermal medium, the charmonium will spontaneously dissociate if the energy e 
obtained from the Schrodinger equation relative to the mass of two open charm mesons 
is positive. The energy e varies with temperature. The critical temperature for the 
spontaneous dissociation (at which e is zero) was determined by Wong [26] to be T d /T c = 
0.99, 0.90, and 0.91 for J/ip, Xci, and Xc2, respectively. Since T/T c < T therm /T c = 0.91 
for any centrality, we can neglect the spontaneous dissociation of J/ip, x&, and Xc2 in the 
pion matter. 

We assume that all charmonia concerned have the same absorption cross section in 
collision with nucleons, <Jnj/iI> = CiVxd = < 7 Nxc2 = a N^' = ^Nabs- Since these cross sections 
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have not been fixed by p— A reactions, we take a value of 4.2 mb. The remaining adjustable 
parameter is r t herm(& = 0) for the most central Pb+Pb collision. The thermalization time 
has a centrality dependence, <5r thcrm ~ l/(a nn n n (b, Tf orm )), where a nn is the nn elastic 
scattering cross section at zero temperature. Since the pion number density is given by 
Eq. (12), the thermalization time after the formation of hadronic matter can be obtained 
by 

Sr theim (b) = Sr thcim (b = o) n ^ b = °> T ^ . (22) 

n n {b, r form ) 

The thermalization time measured relative to the moment of initial collision is then 
^thcrm (b) = r form + 5T thcrm (b). 

The value of the adjustable parameter r thcrm (b = 0) and the transverse velocity 
Kx_i_ are obtained by fitting the latest data of B(J /ip)a(J/ip) /a(DY) from the NA50 
Collaboration [47]. As shown in Fig. 3. the latest data in Ref. [47] are slightly different 
from the earlier data [2,3] in the peripheral and central collisions. Other quantities such 
as the pion number density at thermalization n therm , proper time Tf z , and temperature 
T fz at freeze-out are then calculated and also tabulated in Table I for the central collision 
at b = fm. The quantity n t hcrm = n n (b = 0, r t h C rm) is calculated with Eq. (12). The 
freeze-out time Tf z depends on the centrality. The freeze-out temperature Tf z , determined 
in the most central collision, is the same for hadronic matter yielded at any centrality. 
The temperature between T therm and T fz approximately obey the relation given by Bjorken 
[48]. In the survival probability of total J/^, we include contributions at the level of 58% 
for direct J/ip, 20% for x *, 14.5% for X c2, and 7.5% for f [41]. 

From central to peripheral collisions, pion number density decreases and thermal- 
ization time increases according to Eq. (22). At b = 8.9 fm the pion number density at 
Ttherm is equal to the freeze-out number density and no thermalization is attained. For 
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b > 8.9 fm, the zero-temperature charmonium dissociation cross sections are relevant for 
J/tp suppression. Since the peak dissociation cross sections for J /if), Xci, and Xc2 in colli- 
sion with 7r at T = are about 1.1 mb, 1.6 mb, and 1.9 mb, respectively, the suppression 
of J/ip due to the collision with it is very small. Therefore, the absorption by N — J/ip 
collisions completely dominates J/ip suppression at b > 8.9 fm. For b < 8.9 fm, ther- 
malized hadronic matter sets in with T therm = 0.16 GeV and freezes out at T fz = 0.139 
GeV. The dissociation cross section jumps from 1.1 mb (1.6 mb, 1.9 mb) at T = to 6.5 
mb (4.1 mb, 4.8 mb) at T fz = 0.139 GeV and even higher at T thcrm = 0.16 GeV for J/ip 
(Xci, Xci)- But this sudden rise does not induce a sudden fall of B(J/ip)a(J/ip)/a(DY) 
at Et ~ 40 GeV since the time in which J / ip interacts with pion matter in thermal equi- 
librium increases gradually as the impact parameter decreases from b = 8.9 fm to zero. 
In these calculations, the suppression of ip' has been taken into account (see the following 
paragraph), as ip' comprises about 7.5% of the unsuppressed J /if) yield. 

To study ip' suppression, we note from Ref. [26] that the spontaneous dissociation 
of ip' begins with T d /T c = 0.50. The pion matter freezes out at Tf z /T c ps 0.79 for any 
centrality. Thus ip' spontaneous dissociation takes place throughout the pion matter. For 
a full description of the ip' spontaneous dissociation, a suitable treatment will be needed 
in future work. The effect of ip' dissociation, alternatively, can be accounted for by an 
approximate effective treatment where we assume a constant ip' dissociation cross section 
independent of energy, o^> = 12 mb. As shown in Fig. 4, the experimental data for 
B(ip')a(ip') / 'a(DY) reported in Ref. [49] and exhibited in Ref. [42] can be represented 
by such a 7r — ip' cross section of 12 mb. We shall include this constant tp' dissociation 
cross section in our analysis of the tp' absorption in our model. It can be considered as an 
effective representation of the experimental data to include the effect of ip' suppression. 
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In addition to the spontaneous dissociation, J/ip suppression may receive contri- 
butions from the thermalization of charmonium in hadronic matter. The thermalization 
process leads to a thermal distribution of charmonium, and the states above the dissociate 
threshold can dissociate spontaneously [26, 30]. The thermalization of charmonia depends 
on the dynamic processes such as scatterings of lower charmonium states in collision with 
pions leading to higher charmonium states. The work of Fujii and Kharzeev [31] gave 
the cross sections for ix + J/ip — > ir + ip' less than 0.01 mb in the region of interest. We 
assume implicit generalization of this cross section to other pion-charmonium scatterings. 
Then the thermalization time of charmonia is very long and the charmonia cannot be in 
thermal equilibrium before pionic matter freezes out. Therefore the dissociation by ther- 
malization for J/ip, Xci, and Xc2 have not been included. If we calculate cross sections 
for processes such as n + J /ip — > ix + Xc j in the quark-interchange mechanism, we need 
to interchange the quarks twice and the quark-interchange processes are of higher order. 
The calculations for these cross sections are complicated and are left for a future work. 

To illustrate the role of charmonium dissociations in collisions with pions in hadronic 
matter, we set o"Nabs=0 and calculate B(J/ip)a(J/ip) / 'a(DY) with different assumptions 
on the n — (cc)jls cross sections for the dissociation of the charmonium (cc)jls- In 
Fig. 5 we show the results with n — (cc)jls dissociation cross sections as given by: (1) 
a(n — (cc)jls) for T = [25] (dashed curve); (2) a(ir — (cc) jls)=1 mb (dot-dashed curve); 
(3) o-(tt - {cc) JLS )=2 mb (dotted curve); (4) a(n - (cc) JLS ) of Figs. 1 and 2 for T > 
(solid curve). The flat curve below E T = 28 GeV comes mainly from the collisions at 
b >8.9 fm where no thermalization is attained. This flatness arises as the 7r-charmonium 
dissociation cross section is small and the pion number density is not high at b > 8.9 
fm. Indeed, the dashed curve gives a small J/ip suppression if only the 7r-charmonium 
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dissociation cross sections at T = are used. Even the J/ ip suppression obtained by using 
1 mb cross section is larger than that for the cross sections at T = 0. It is obvious that the 
suppression obtained from the cross sections at higher temperatures is greater than that 
from using a 2 mb cross section. This is related to the larger charmonium dissociation 
cross sections at higher temperatures as shown below. 

The cross sections in Figs. 1 and 2 depend on the temperature and ix — (cc) JLS 
relative momenta. What are the average cross sections? Fig. 6 exhibits cross sections 
averaged over the pion and (cc)jls momenta in the collision at b = fm, which is defined 



as 



E v (Pp / (2tt) 3J " v ' I 7 E y d 3 p 'J (2vr) 3 ' 
The average cross sections (cr n j/^), (cr nXcl ), and (o- nXc2 ) decrease with time. At the ther- 
malization time, they are 3.47 mb, 2.74 mb, and 2.95 mb, respectively. At 2.9 fm< r <4.25 
fm, (a w j/^) is larger than (cr^xd)- ^ 2-9 fm < r <3.4 fm, (cr n j/^) even surpasses {o- nXc2 ). 
These are not surprising since the peak value of J/ip dissociation cross section exceeds 
that of Xci (Xcs) f° r T/T c > 0.55(0.65) and the temperature of hadronic matter before 
freeze-out is T/T c > 0.79. 

We further calculate the quantity 



W= dr{v Tel cr)/(Tf z - r thcrm ). (24) 



The pr and xf dependences of v re \a are shown in Figs. 7 and 8. These results indicate 
that the average absorption effects due to J/ip, Xa, and Xc2 are not significantly different. 
The above discussion will be useful to understand data from the NA60 Collaboration 
where the Xcj suppressions will be measured. 

C. J/ip Suppression at Large E T 

17 



In addition to the general suppression over a broad region of transverse energy E T , 

the data for J ftp suppression in the large transverse energy region with E T > 100 GeV 

are also of special interest. For E T > 100 GeV, we consider the effects arising from the 

large Et fluctuation. Since D(b, Et) has not been constrained to give Et by the average 

jd%E T (b)D(b,E T ) 
{Et) = J d 2 bD(b, Et) ' (25) 

the measured E T is larger than (Et) for E T > 100 GeV, i.e. the so-called large-E^ 
fluctuation happens [15,16]. In our calculations, the fluctuation increases not only the 
pion number density by the replacement n n (b, r) — > n n (b, t)Et/ (Et) but also the lifetime 
of pion matter by V(b, Tf z (fe)) = V( y b,T t i ierm (b))n 7T (b,T t \ 1CTm (b))/n{ 7i . The two changes lead 
to the solid curve in Fig. 3 decreasing continually beyond E T = 100 GeV. 

The large E T fluctuation does not increase the initial temperature since we do not 
consider here the QCD phase transition to a QGP. The spontaneous dissociations of J ftp, 
Xd, and Xc2 can still be neglected. The factor E T /(E T ) shows the fluctuation increases 
smoothly from E T = 100 GeV. Therefore we see that the solid curve does not drop 
rapidly in the large- Et region and does not pass the rightest experimental point within 
the error bar. This behavior is similar to the theoretical J ftp suppression obtained by 
Capella et al. [15]. If experimental data at Et > 120 GeV can be provided with good 
statistics by additional NA60 measurements [50] , then we need compare our result to the 
experimental data in this region to see the possibility of additional suppression due to the 
QGP. 

VI. CONCLUSIONS 



We have studied 7r-charmonium dissociation cross sections and the J ftp suppression 
in Pb+Pb collisions at the CERN-SPS. We use dissociation cross sections calculated in 
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the Barnes-Swanson quark-interchange model [32,33] with the potential obtained from the 
lattice gauge results. Making reasonable assumptions about the absorption scenario, we 
find that the J/ip production cross section calculated with these dissociation cross sections 
can describe the general features of the anomalous J/ip suppression data. An important 
element of the agreement arises from the increase of the dissociation cross sections as 
the temperature increases. For the large E T region at E T > 100 GeV, the anomalous 
additional suppression can be explained by the fluctuation of E T or multiplicity. 

The 7r-charmonium dissociation cross sections depend on the temperature and the 
kinetic energy, but the averages v re \a in hadronic matter for ir + J/tp, n + Xci, and tt + Xc2 
dissociation cross sections are nearly the same. Thus, the dissociations of direct J/ip, Xci, 
and Xc2 give nearly the same contributions to the suppression of the measured J/ip. 

We have chosen the constant absorption cross sections of charmonia on nucleons. 
However, the absorption cross sections are expected to depend on the CM. energy of 
charmonium and nucleon. Determining the nucleon absorption cross sections is an im- 
portant task. It will be of interest to evaluate the N — J/ip dissociation cross sections in 
the Barnes-Swanson quark-interchange model [32,33]. 
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Table 1: Values of <T Na bs, V ex ±, r t herm, n t hcrm, r fz , and T fz . 



o- Nabs (mb) 


Kx±(c) 


Ttherm(Wc) 


n t hcrm(fm d ) 


r fz (fm/c) 


T fz (GeV) 


4.2 


0.75 j 


2.9 


0.88 


4.44 


0.139 




Figure 1: J/ip dissociation cross sections for various temperatures as a function of Eke- 
From left to right are T/T c =0.45, 0.35, 0.25, 0.15, 0.05 for these curves without labels. 
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10 




Figure 2: Upper and lower panels show Xci an d Xc2 dissociation cross sections for various 
temperatures as a function of E KE . From left to right are T/T c =0.35, 0.25, 0.15, 0.05 for 
these curves without labels. 
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Figure 3: The experimental data are from Ref. [47]. The solid curve is the theoretical 
result while UNabs = 4.2 mb. The dashed curve is obtained with no large-Er fluctuation 
considered. 
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Figure 4: Solid curve obtained with 0N a bs = 4.2 mb is compared with the 1996 data of 
a{ip')/a(DY) ratio. 
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Figure 5: The same experimental data as Fig. 3. All curves have a Nabs = mb. The 
solid, dashed, dotted, and dot-dashed curves are obtained by using the cross sections at 
T > and T = 0, the two constant cross sections cr n j/^ = cr^xcj = = 2 mD an d 1 mb, 
respectively. 
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x (fm/c) 

Figure 6: The solid, dashed, and dotted curves are (cr n j/^), {cr WXcl ), and (ov Xc2 ), respec- 
tively. 
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Figure 7: The solid, dashed, and dotted curves are v ie ia n j/^, freicVxci , an d v lc \aT, Xc2 at 
xf = 0, respectively. 
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